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Abstract 

An extended study is performed of geometrical and kinematical assumptions 
used in calculations of the neutrino oscillation phase. The almost universally em¬ 
ployed ‘equal velocity’ assumption, in which all neutrino mass eigenstates are pro¬ 
duced at the same time, is shown to underestimate, by a factor of two, the neutrino 
propagation contribution to the phase. Taking properly into account, in a covari¬ 
ant path amplitude calculation, the incoherent nature of neutrino production as 
predicted by the Standard Model, results in an important source propagator con¬ 
tribution to the phase. It is argued that the commonly discussed Gaussian ‘wave 
packets’ have no basis within quantum mechanics and are the result of a confused 
amalgam of quantum and classical wave concepts. 
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1 Introduction 


The first published calculation of the phase of neutrino oscillations, on the assump¬ 
tion that neutrinos are massive, but ultra-relativistic, was made in 1969 by Gribov and 
Pontecorvo [jTj. No details of the method used to obtain the phase were given. Later 
publications used this result for phenomenology. In 1976 a paper was published [|| in 
which an oscillation phase a factor of two smaller than the Gribov-Pontecorvo result was 
obtained. No comment was made, then or later, on this discrepancy, which remained 
unnoticed before a paper written recently by the present author [|J in which more biblio¬ 
graphic details of this ‘missing controversy’ can be found. Since the first derivation of the 
‘standard’ oscillation phase, as the result of Reference [|J will be referred to in this paper, 
of the order of a hundred papers have been published devoted to the quantum mechanics 
(QM) of neutrino oscillations. 

The present paper has three aims. The first is to carefully review the kinematical 
and geometrical assumptions of previous calculations, particularly with regard to their 
consistency in terms of an expansion in powers of the neutrino masses. The second is 
to present a calculation of the oscillation phase, for sources at rest, using a covariant 
Feynman path amplitude formulation f|, [| of QM 0 The third is to critically discuss 
treatments, in the literature, of the QM of neutrino oscillations, in the light of both 
the kinematical and geometrical findings of the present paper and the physical insights 
provided by the path amplitude treatment. 

The calculations presented here are only a sub-set of those to be found in a previous 
paper by the present author ||, results of which are also summarised in a short letter j6j. 
The calculations in Reference |3j include neutrino oscillations following decay at rest 
pions, muons and /3-unstable nuclei, muon oscillations following pion decay at rest, and 
neutrino and muon oscillations following decay in flight of ultra-relativistic pions. Also 
discussed in Reference Q are different mechanisms (both coherent and incoherent) that 
may contribute to damping of the oscillations. An attempt has been made, in the present 
paper, to give a more pedagogical presentation in which the crucial underlying physical 
assumptions used are spelt out as clearly as possible. 

Reference || contains already a extensive critical review of the previous literature. In 
the present paper, the further discussion has been largely motivated and orientated by the 
reports of several anonymous reviewers of References |3], ||]. Since writing these papers I 
became aware of the important work of Shrock [[7. §] who pointed out, more than twenty 
years ago, the incoherently nature of the weak decay processes in which neutrino mass 
eigenstates are produced. This was indeed one of the crucial assumptions for the calcula- 

-Reyman’s original work on this subject cited in References ||, || concerned only non-relativistic QM, 
but the fundamental formula (5.1), on which the work presented in the present paper is based, holds also 
in relativistic QM. The approach should be more properly termed a ‘path amplitude’ rather than a ‘path 
integral’ one as no attempt is made to evaluate the latter in closed mathemetical form, as is done, for 
example, in the discussion of bound state problems. 

2 The word ‘incoherent’ here means that the different mass eigenstates are produced in independent 
physical processes, not in a coherent ‘flavour eigenstate’ that is a superposition of the mass eigenstates. 
A ‘coherent source’ is defined as one in which the amplitudes corresponding to different decay times have 
a well-defined phase difference. 
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Tabic 1: Values of the flavour-mass mixing amplitude ( p\q) 


tions presented in References ||, ||. I also now fully understand the (simple) reason why 
the standard oscillation phase has been obtained in essentially all published derivations 
since Reference [|| . Indeed, as will be demonstrated, the wrong (standard) instead of the 
correct (Gribov-Pontecorvo) neutrino propagation phase has been universally obtained 
due to the universal failure to use the correct description of the weak neutrino production 
processes as formulated by Shrock [7:. |5j . 

For simplicity, only two-flavour mixing will be considered, so that the mixing ampli¬ 
tudes: {i\a) (where i (= 1 , 2 ) refers to neutrino mass eigenstates zq,z/ 2 and a (= e,/i) 
refers to the flavour of the charged lepton participating in the weak interaction process) 
may be taken to be real. Unitarity then enables all mixing amplitudes to be described in 
terms of a single angular parameter, 6 , as shown in Table 1. Note that there are no ‘neu¬ 
trino flavour eigenstates’. The amplitude (1 |/x), for example, is the amplitude to produce 
the neutrino mass eigenstate zq in association with (or in the decay of) a muon. Thus 
the flavour-mass mixing amplitudes fix the strength of the charged current couplings to 
the neutrino mass eigenstates 0, |[. For example, in the decays: 7 r + , K + , B + —> jU+zq 
the leptonic charged current is: Ff 7 P (l — 75 )f/ p ?./i (i = 1 , 2 ) to be compared with the 
quark currents: q^ p (l — 75 )V uq u ( q = d,s,b ), where U fn = U itl = (l|/u) is an element of 
the Maki-Nakagawa-Sakata (MNS) |9[| lepton flavour/mass mixing matrix, and V uq is an 
element of the Cabibbo-Kobayashi-Maskawa (CKM) JlOj flavour/mass mixing matrix in 
the quark sector. The couplings of antineutrinos to charged anti-leptons are the same as 
those of neutrinos to charged leptons in the two flavour mixing case. 


The amplitude to produce a charged lepton of flavour (3 by the interaction of a neutrino 
mass eigenstate i produced in a decay involving a charged lepton! of flavour a may, in 
general, be written as: 

Aflfl <- a) = A o exp(-iA0 i )(/3|i)(i|a') (1.1) 


The probability to observe a charged lepton of flavour (3, in the case that the neutrino 
mass eigenstates are not distinguished, is then, by superposition: 


P/3a — \Ai(/3 <— a) + A 2 {(3 <— a) | 2 

= | A 0 | 2 [((/?| 1 )( 1 |«» 2 + (</?| 2 )< 2|«)) 2 

+2(/3|l)(/?|2)(l|a)(2|a) cos(A0i - A0 2 )] 


( 1 . 2 ) 


The ‘neutrino oscillation phase’ 0 i2 is defined as 


012 = A0i — A0 2 (1-3) 

3 In the decay process, this charged lepton, together with the neutrino m;, form a charged current that 
couples to a real or virtual W boson. 
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The Gribov-Pomeranchuk value of the oscillation phase is: 


whereas the standard result is: 


rp A m 2 L 
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(1.4) 
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(1.5) 


In Eqns(1.4) and (1.5), Am 2 = m\ — m| where m\ and m 2 are the neutrino masses, L 
is the source-detector distance and and p v is the measured neutrino momentum. Units 
with h — c — 1 are used throughout this paper. The word ‘neutrino’, without further 
qualification, stands for ‘neutrino mass eigenstate’. The plan of the paper is as follows: In 
Sections 2, 3 and 4 various kinematical and geometrical approximations that have been 
used in neutrino oscillation calculations are discussed. In particular it will be examined 
whether particular approximations retain, or not, all the leading order 0(m 2 ) terms, 
or contribute only negligible 0(m 4 ) corrections. The Lorentz invariant nature of the 
oscillation phase is fully exploited in these purely mathematical considerations. Section 
5 describes the calculation of the oscillation phase using the covariant Feynman path 
amplitude method. Section 6 is devoted to a discussion of previous treatments of the QM 
of neutrino oscillations. In particular, the physical basis of widely used Gaussian wave 
packet models is questioned. Section 7 contains a brief summary and outlook. 


2 Lorentz Invariant Plane Wave Propagation 


In this case, the one-dimensional propagation amplitude for a particle of mass m, : over 
space and time intervals: A Xi and A ti is given by: 


P(Axi,Ati,mi) = P 0 exp{-i[EiAti - p-A'Xi}} 

= P 0 exp {—zm,;Arj} 

= Pq exp {—iAcfti} (2.1) 

Here At* is the proper time interval, in the rest frame of the particle, i, corresponding to 
space and time intervals Axi and At, in the laboratory frame, so that: 

At 2 = (A Uf - (A Xi f (2.2) 

The relativistic relation between energy, momentum and mass: 




(2.3) 


implies a group velocity vf for the travelling wave represented by Eqn(2.1) given by: 


v? = 


dEj 

dpi 


Pi_ 

Ei 


(2.4) 


In the following, it is assumed that the neutrino mass eigenstates whose propagation is 
described by Eqn(2.1), have velocity Vi = vf. As correctly emphasised in Reference [[TI 
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particle oscillation experiments actually measure the spatial, not the temporal, depen¬ 
dence of quantum mechanical interference effects. Therefore, throughout the kinematical 
discussions in the present and following Sections, a fixed distance, A ay = L, is assumed 
between the source particle (at rest) and the detection event. With the additional as¬ 
sumption that the propagating neutrino mass eigenstates are on-shell particles, several 
different exact expressions may be written for the Lorentz-invariant phase increment: 


A fa = rriiATi 


rriiAti 


7» 


rrtfAti rnfL 

Ei Pi 


(2.5) 


where, in the last member of Eqn(2.5), the relation: A ti = L/vi has been used, and 
7 = l/\/l — v 2 = E/m is the usual parameter of Special Relativity. 


The value of the oscillation phase: = A<p\ — A0 2 obtained using different kinemat¬ 

ical approximations will now be considered. For this it is useful to introduce the quantity 
Po which is the limiting value of pi or Ei as m,; —> 0. The following expressions for p t and 
Ei may then be written: 


Pi 


Ei 



Po {1 + 



1 + R 


(1 ~R 2 m) 2 

1 

1 - B 2 

J. ± L m 


0(m 4 ) 


( 2 . 6 ) 

(2.7) 


where 

Po = ^(l-C) ( 2 . 8 ) 

Here ms is the mass of the decaying source particle and R m = ms/ms where rrin is the 
mass of the particle (or system of particles) recoiling against the neutrino in the decay 
process. It is also convenient to introduce the average kinematical quantities for the 
neutrino mass eigenstates: 


m 

P 

E 


(mi + m 2 ) 
2 

(Pi + Pr) 

2 

Jp 2 + m 2 


(2.9) 

( 2 . 10 ) 

( 2 . 11 ) 


If, as is usually done in the literature, common values of both Ax and At are assumed 
for both mass eigenstates, then: 


Vi = v 2 = v = 


Ax 

A t 


P_ 

E 


and 


7i = 72 = 7 = 


At _ E 
At 


m 


( 2 . 12 ) 


(2.13) 


It must be remarked at once that such an ‘equal velocity’ hypothesis is in contradiction 
with the fundamental definitions in Eqns(2.1)-(2.4). For logical consistency, if Ax and 
At are assumed to be equal for the two mass eigenstates, then also Ei and pt in Eqn(2.1) 
should be replaced by E and p. Then, evidently, no neutrino oscillations would be possible. 
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In fact the equal velocity hypothesis assumes not only Eqn(2.12) but also that the energy 
and momentum appearing in the 4-vector product in (2.1) are in accordance with energy 
momentum conservation in the decay process. Since vf = p,JE t p v, there is then an 
internal contradiction. Also, because of Eqn(2.12), the corresponding space and time 
intervals do not respect the equations: vf = Axi/Ati, i.e. the neutrinos do not move 
along classical trajectories specified by the group velocities vf. Thus the equal velocity 
hypothesis respects energy-momentum conservation, but not the constraints of space- 
time geometry, and contains an internal contradiction. The formula (2.1) combines the 
quantities (Aay, At*) from space-time geometry and (E l} p,) from particle kinematics. 
The equal velocity hypothesis then treats, in an inconsistent manner, these two sets 
of quantities. Introducing the notation A (f>i(v) for the approximate phase increment 
calculated with fixed velocity v, according to Eqns(2.9)-(2.13), then: 


A (j>i(v) = mjAr 


rriiAt 

7 


m 2 At 

E 



(2.14) 


Now denoting by the superscripts I, II, III and IV the values of 0 12 (f) calculated using 
the successive members of Eqn(2.14) the following results are obtained: 
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(2.15) 

(2.16) 

(2.17) 

(2.18) 


The phases and are found to be equal and a factor of two larger than either 0{ 2 or 
0 ( 2 - The different results obtained indicate that some inconsistent approximations must 
have been made in the passage from the exact expressions for the phase increment in (2.5) 
to the approximate ones in (2.14). To see just where the inconsistency has occured, and 
which results, / and II or III and IV, are correct, the calculation is now done using the 
exact formulae of Eqn(2.5) and, for ease of comparison with the approximate expressions 
in Eqn(2.14) the neutrino masses are written as: 


m, = m + Si 


(2.19) 


where 


A — —A 2 



( 2 . 20 ) 


Thus Si are correction terms, relating the approximate relations in Eqns(2.14) to the exact 
ones in Eqn(2.5). Denoting by 0i2 (exact) the interference phase calculated using Eqn(2.5) 
then: 


0 i 2 (exact) 
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( 2 . 21 ) 
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Thus the result obtained is in agreement, at 0(m 2 ), with the approximate calculations in 
cases III and IV. Using now Eqns(2.19) and (2.20): 


0 i 2 (exact) = 


(m + <5i) 2 

(m + 5 2 ) 2 

J — 

(m + <5i) 2 - (m + 5 2 ) 2 

Pi 

P2 

Lj — 

P 


K-m|) L + (mi - m 2 )m L + < 

2 p p 


L + O/m 4 ) 

( 2 . 22 ) 


Comparing the second term on the RHS this equation with Eqns(2.15) or (2.16) it can 
be seen that the correction terms in Eqn(2.19) have the effect of doubling the approx¬ 
imate results 0i 2 (u) and 4>[l(v). Thus important O(m 2 ) contributions are neglected in 
Eqns(2.15) and (2.16). The reason why the constant velocity approximation is so poor in 
Eqns(2.15) and (2.16) becomes evident on inspection of Eqn(2.13). Unlike in the case of 
neutral kaons or b-mesons, the neutrino masses may be widely different, so that although, 
for highly relativistic neutrinos, the velocities of the mass eigenstates may be very similar, 
this is not the case for the Lorentz-y factors in Eqn(2.13) that are inversely proportional 
to the neutrino masses. The important difference between the exact and approximate cal¬ 
culations occurs in the very first member of Eqn(2.15). Indeed, as is shown by Eqns(2.21) 
and (2.15): 


0i2 (exact) — 0{ 2 ( y ) = m-iAri — m 2 Ar 2 — (mi — m 2 ) At 


{m\ — ml) 

2 p 


L + 0(m 4 ) (2.23) 


In fact, in Eqns(2.17) and (2.18), the approximate Lorentz-y factor y = E/m is replaced 
by the still approximate but mass dependent y-factors E/mi which differ from the exact 
ones Ei/mt by corrections of only 0(m 2 ). Thus, up to corrections of O(m 4 ), the same 
result is obtained as the exact formula. On the other hand the difference between 1/y, 
used in Eqn(2.16), and 1/yi or l/y 2 used in Eqn(2.5), is of O(m) and so cannot be 
neglected in an 0(m 2 ) calculation of the oscillation phase. 


It is instructive to calculate the oscillation phase directly in the laboratory system, 
using Eqn(2.1), rather than by Lorentz transformation from the neutrino center of mass 
to the laboratory system as in Eqns(2.5) and (2.14). Considering first the equal velocity 
case: 


0i2 4B ( u ) = ( E i - E 2 )At - (pi - p 2 )L 

= [(E 1 -E 2 )^-(p 1 -p 2 )]L (2.24) 

Using the relations: 

2 

Ei = y/pi + rrif = Pi + + 0{m 4 ) (2.25) 

^Vi 

_2 

777 

v = 1-— 2 + 0 (m 4 ) (2.26) 


and noting that, from Eqns(2.6) and (2.10) that p t = p+0(m 2 ), Eqn(2.24) may be written 
as: 


^ B iy) 
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2 p 
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(2.27) 
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in agreement with Eqns(2.15) and (2.16). The exact formula gives: 


</>i 2 B (exact) = 


P 

in agreement with Eqns(2.17) and (2.18). It can be seen that: 

4>i^ B (exact) — 0i2 4 ' B (t) = (EiAti—E 2 At 2 ) — (Ei—E 2 )At = -——— L+0(m 4 ) (2.29) 

2 p 

Again, the ‘equal velocity’ formula is seen to neglect O(m 2 ) contributions that double the 
value of the oscillation phase. 

The neglect of important O (m 2 ) terms when the equal velocity hypothesis is made 
especially transparent on rewriting Eqn(2.21) in the following way: 


(EiAE - E 2 At 2 ) - (pi - p 2 )L 
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E\ E 2 

- Pi -b P2 

Vi v 2 


m t 




Pi P 2 

(m\ — m|) 


L + 0(m 4 ) 


(2.28) 


0i2 (exact) = m\Ari — m 2 Ar 2 = (mj — m 2 )Ar 1 + m 2 (Ari — A r 2 ) (2.30) 


Since 


Eqn(2.30) can be written as: 


A Ti = 


miL rriiL 3 

-=-h 0(m ) 

Pi Po 


(2.31) 


cj)i 2 (exact) = -— -—[mi + m 2 \L + 0(m 4 ) (2.32) 

Po 

The equal velocity hypothesis now implies: 

(a) Neglect of Ari — Ar 2 i.e. of the second term in the square brackets in Eqn(2.32). 

(b) Making the replacement: rrii/p 0 —> m/p in the remaining term of Eqn(2.32). 


It is now obvious that (a) and (b) neglect (different) 0(m 2 ) terms; a term ~ m 2 (mi — m 2 ) 
for (a) and a term ~ (mi — m)(mi — m 2 ) = (mi — m 2 ) 2 /2 for (b). The sum of these 
neglected terms is ~ (m 2 — m|)/2 = m(mi — m 2 ), equal to the terms retained by the 
equal velocity hypothesis. 

Thus the essentially universally employed ‘equal velocity’ assumption results in the 
neglect of 0(m 2 ) terms in the calculation of the vacuum oscillation phase. The existence 
of these terms explains the factor of two difference between the original calculation of the 
oscillation phase by Gribov and Pontecorvo (which is the correct 0(m 2 ) result) and the 
‘standard phase’ given by all calculations that make the ‘equal velocity’ assumption. 

The calculation of the oscillation phase is now repeated assuming, instead of equal 
velocities for the neutrino mass eigenstates, either equal momenta or equal energies. For 
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equal momenta the following relations hold: 


L 

= A ti(p)vi(p) 

(2.33) 

Ei(p) 

= \Jp 2 + nrif 

(2.34) 

Mp) 

p 

Ei(p) 

(2.35) 

7 i(p) 
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(2.36) 


These yield, for the invariant phase increment: 


A <f>i(p) = rriiATiip) 


mjAti(p) 
7 tip) 


m%Atj(p) 

Ei{p) 



(2.37) 


Comparing with Eqn(2.5) it can be seen that each member of the two equations differ 
only by terms of 0(m 4 ). Thus, at O(m 2 ) the Gribov-Pontecorvo result (2.21) is obtained 
for the oscillation phase. For the equal energy case: 


L 

= A ti(E)vi(E) 

(2.38) 

ME) 

= Je 2 - 

(2.39) 

Me) 

n(E) 

E 

(2.40) 

ME) 

E 

(2.41) 

m 

rriiAME) = 

rriiAti(E) m 2 Ati(E) m 2 L 

ME) E ME) 

(2.42) 


Again, all terms in last equation differ from the corresponding ones in the exact expres¬ 
sion (2.5) only by terms of 0(m 4 ) so that the oscillation phase agrees with the value in 
Eqn(2.21). It is easily seen that this result is also obtained in the case of either equal 
momenta or equal energies if the oscillation phase is calculated directly in the laboratory 
system as in Eqn(2.28). 


It may be shown, in a similar way, that once a constant value of At is assumed, 
the standard oscillation phase is obtained independently of whether the energy and mo¬ 
mentum are calculated on the assumption of exact energy-momentum conservation, or 
whether the equal momentum and different energy or different momentum and equal 
energy hypotheses are used. 


Thus, the ‘standard’ oscillation phase is obtained only in the case of the ‘equal veloc¬ 
ity’ hypothesis. In all cases: ‘equal velocity’, ‘equal momentum’ and ‘equal energy’ an 
unphysical assumption is being made (i.e. one that violates conservation of energy and 
momentum). The study just presented shows that in the latter two cases the error of the 
approximation is only of 0(m 4 ) in the oscillation phase, and so is negligible at 0(m 2 ). 
However, in the former case, the neglected terms are of O(m 2 ) and so lead to an incorrect 
result at this order. 


Finally, in this Section it will be found of interest for the subsequent discussion of 
independent temporal or spatial propagation to calculate the separate temporal (T) and 



spatial (S) parts, in the laboratory system, of the Lorentz-invariant oscillation phase (pn- 
tf 2 AB (T) = (EiAti - E 2 At 2 ) 
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= \Pl - P2 H- 

V p i 

(1 + R 2 m ) 


1717 , 
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1717 , 


Po 


-L + 0(m 4 ) 


(2.43) 


= -(p,-p 2 )L= {1 + R2J (mi-ml)L + 0(m i ) (2,44) 

4p 0 

On adding (2.43) to (2.44) it can be seen that the Lorentz invariant oscillation phase 
originates entirely from the temporal part, the spatial part being exactly cancelled by a 
term in the temporal one. 


3 Gaussian Wave Packet Propagation 


Very often in the literature, following Reference the plane wave propagator of 
Eqn(2.1) has been modified by the introduction of a multiplicative Gaussian spatial wave 
packet: G(Ax, At, Vi, a x ): 

P WP (Ax, At, mi, Vi, a x ) = G(Ax, At, v h a x )P(Ax, At, (3.1) 

where P(Ax, At, m;) is given in Eqn(2.1) and : 

'Ax — UjAO 21 


G(Ax, At,Vi,a x ) = (V2na x ) 2 


exp 


2(T A 


(3.2) 


In this formulation the intervals Ax and At are assumed to be the same for all mass 
eigenstates, so that the equal velocity hypothesis is tacitly introduced. To calculate the 
oscillation phase, 0^2 P an integration over the transit time in the range —00 < At < 00 
is performed, leading to the result fl^| : 

— - (Pi - P 2 ) 


% P 


E\ 


where 


vi + Un ^ 1 / m ? m 2 \ 4 , 

v = — -- = l- - -^ + -£ + O m 4N 

v 1 + v 2 2 [pj ~ 2 1 y 


pi 


(3.3) 


(3.4) 


It may be noted that the effective average velocity v is lower than that of either mass 
eigenstate. This is presumably due to the contribution of the unphysicali negative region 
in the integration over At. Using the relation: 

,2 


Ei = Pi 


2pi 


+ 0(m 4 


(3.5) 


4 Note that the Gaussian variation of At, for fixed source detector distance At = L , assumed in this 
calculation, is at variance with the known space-time structure of the decay and detection events, i.e. an 
exponential distribution of decay times followed by a mass dependent propagation time At,;. 
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Eqn(3.3) may be written: 




(Pi - P 2 ) + 


1 mi mi \ 

2 + 



(pi -P 2 ) f + 0(m 4 ) 


L + 0(m 4 ) 


(m^ — m^) 


L + 0(m 4 ) 


(3.6) 


where Eqn(2.6) has been used to write Pi in terms of p 0 and mThe ‘wave-packeted’ 
propagator thus gives the same oscillation phase at 0 (m 2 ) as the unmodihed invariant 
plane wave propagator when the same equal velocity hypothesis is made. In both cases 
important 0 (m 2 ) terms are neglected in comparison with the exact calculation. 


In a more elaborate recent calculation [|K|] where wave packets are associated not only 
with the neutrinos, but also with all particles participating in the production and decay 
processes, a formula for the phase increment identical to (3.1) was obtained except that 
the width of the spatial wavepacket a x is replaced by a parameter ■// that depends on 
the widths of the wavepackets of all participating particles. Since the oscillation phase 
(3.3) does not depend on this parameter the same ‘standard’ result is obtained as in 
Reference |llj. The result obtained with different mean velocities, v given by Eqns(2.9)- 
(2.12) and v given by Eqn(3.4) are the same. Indeed, it is clear from the derivation of 
(3.6) from (3.3) that a mean velocity (v) defined by any formula of the type: 


{v) 


am\ + /3mim 2 + qm| 4 

1- 35 -b Oim ) 

p 


(3.7) 


where a, (3 and 7 are arbitary coefficients of order unity, will yield the result (3.6). 
However, in the derivation of the same formula in Eqns(2.15) ,(2.16) this result requires 

that m = (mi + m 2 )/2. If instead the value m = \Jm\ + mf, suggested by Eqn(3.4), is 
used, then a different result is obtained: 



(? 7 ii - m 2 )(\Jm\ + ml)L 


p 


+ 0(m 4 


(3.8) 


Thus, unlike in the case of the invariant plane wave, the wave-packet treatment gives 
inconsistent results for the oscillation phase in the neutrino centre-of-mass and laboratory 
systems. 


4 Independent Temporal or Spatial Propagation 


Many descriptions of neutrino oscillations (especially concise presentations in review 
articles concerned mainly with the description of experimental results, or theoretical mod¬ 
els of neutrino masses) do not use a Lorentz invariant description but consider instead 
independent temporal (TE) or spatial (SE) evolution of the wavefunctions of the neutrino 
mass eigenstates. It is in this way that the standard formula for the oscillation phase of 
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Eqn(1.5) was first derived in Reference ||. A simple Schrodinger time evolution is used 
giving the temporal propagator: 


P T E(At,mi) = Pq exp [-iEiAt] 


(4.1) 


or 

A <j>i(TE) = EiAt (4.2) 

To relate the time interval At (assumed to be the same for both neutrino mass eigenstates) 
to the experimentally measured source-detector separation L, it is further assumed that 
the neutrinos are ultra-relativistic and so move at the speed of light, i.e. At = L. To 
derive the oscillation phase the eigenstates are then assumed to have equal momentum, 
p v , but different energies. As first pointed out in Reference Jl4| this assumption is in con¬ 
tradiction with energy-momentum conservation in the decay process. This is easily seen 
by inspection of Eqns(2.6) and (2.7) above. Indeed, the violation of energy-momentum 
conservation is an 0(m 2 ) effect. Still, making this assumption: 

^=P, + ^ + 0(m 4 ) (4.3) 

so that Eqn(4.2) gives: 

MTE) = AMTE) - A<j> 2 (TE) = (m? ~ ^ + 0(m 4 ) (4.4) 

2 p v 

This calculation is wrong, firstly because it does not respect Lorentz invariance and sec¬ 
ondly because the kinematical approximations made (neutrinos of equal momenta moving 
at the speed of light) neglect important 0(m 2 ) terms, as can be seen by comparing the 
exact expression of the temporal part (f>i£ B {T) in the laboratory system of the Lorentz 
invariant oscillation phase in Eqn(2.43) with Eqn(4.4). 


Returning now to the question of the compatibility of Eqn(4.1) with special relativity, 
such a formula is a reasonably good approximation to the exact, Lorentz-invariant, result 
in the non-relativistic (NR) limit where pi -C m*. In this case: 

A tf AB {NR) = EiAU = ^(l + £) (4.5) 

Pi V m iJ 

so that, up to a fractional correction p 2 /m|, the non-relativistic result is the same as 
the exact one (2.5), and the same result, Eqn(2.21), is obtained for <^ 12 . The different 
value obtained in Eqn(4.4) evidently is a consequence of the ‘equal time’ assumption 
At i = A t 2 = At, which results in the neglect of important 0(?n 2 ) terms. It should be 
stressed that the use of Eqns(4.4) or (4.5) to describe ultra-relativistic neutrinos is wrong 
in just the same way that the use of the non-relativistic formula for the total energy of a 
particle, E = m + T, where T = mv 2 / 2, is wrong when v ~ 1. 


Other authors [O, ITR ITR] motivated by the correct observation that flavour oscillation 


experiments actually measure quantum interference effects as a function of space, not time, 
have proposed to describe neutrino oscillations in terms of only a spatial evolution of the 
wavefunction: 


Pse = P 0 exp[ipiL\ 


(4.6) 
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or 


A faiSE) = - Pi L (4.7) 

Now, making the additional assumption (incompatible with energy-momentum conserva¬ 
tion in the decay process) of equal energies, E v , so that: 


Pi E v 


2 E v 


+ 0(m 4 ) 


(4.8) 


Eqn(4.7) gives: 


(pn(SE) = AMSE)-A<j) 2 (SE) = 


(ml — ml) L 

2 K 


+ 0(m 4 


m 


m 2 2 )L 


2 p u 


+ 0(m 4 ) 


(4.9) 


in agreement with Eqn(4.4). As in the case of independent time evolution, Eqn(4.7) is 
wrong, firstly, because it does not repect Lorentz invariance, and secondly kinematical 
assumptions are made (equal energies and different momenta) leading to a result differ¬ 
ing by 0(m 2 ) terms from the spatial part of the exact invariant result, as can be seen 
by comparing Eqns(4.9) and (2.44). Unlike for the case of time evolution, there is no 
kinematical domain in which an equation similar to Eqn(4.7) is a good approximation. 
It is tantamount to performing a kinematical calculation with 4-vectors, the temporal 
components of which have all been arbitarily set to zero. In all cases it yields a very bad 
approximation. Actually, to derive Eqn(4.6) [|ll], [15], [[(J the Lorentz invariant plane wave 
was written down and equal velocities, i.e. A t\ = A t 2 = At, as well as equal energies 
were assumed for both mass eigenstates so that the temporal contributions cancel exactly 
in the interference phase. This is equivalent to assuming the purely spatial propagation 
of Eqn(4.6). The physical arguments given in References [|lT], |U| [Rj to justify the equal 
energy hypothesis are further critically examined in Section 6 below. 


It should be noted that both the TE and SE hypotheses contain a further internal 
contradiction, different from that of the equal velocity hypothesis discussed in Sections 3 
and 4 above. In the latter, the correct values of E t and pi, as given by energy-momentum 
conservation, are assigned in the plane wave propagator, only the geometrically inconsis¬ 
tent values Afi = At 2 = At are assumed. For the TE and SE hypotheses, not only are the 
values of AU inconsistent with the relation A ti = L/vi = LEi/pi, but energy-momentum 
conservation is also violated by the ‘equal momentum’ or ‘equal energy’ hypotheses. The 
latter assumptions, however, only generate shifts of 0(m 4 ) in the oscillation phase as 
compared to the result obtained assuming exact energy-momentum conservation. 


5 The Quantum Mechanics of Flavour Oscillations 


In this Section, a general discussion of the underlying principles of flavour oscillation 
calculations in the Feynman path amplitude formulation of QM is first given before deriv¬ 
ing the oscillation phase for a simple two-flavour problem with the source particle at rest. 
More refined calculations performed in a similar manner may be found in Reference || . 
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After the pioneering conceptual work of Planck, Einstein and Bohr, QM was for¬ 
mulated in several independent ways; most importantly, this was done by Heisenberg, 
Schrodinger and Feynman. The Heisenberg (Matrix Mechanics) and Schrodinger (Wave 
Mechanics) formulations originally competed for the description of atomic physics. How¬ 
ever, the greater similarity with classical physics, the less abstract nature and superior 
calculational power of the Wave Mechanics approach resulted in this soon being adopted 
as the standard one, both in atomic physics research and in text books on non-relativistic 
QM. Wave Mechanics is particularly well adapted to subjects such as atomic, molecular 
and nuclear physics, where the solution of bound state problems plays an important role; 
a domain that may be called ‘Quantum Statics’. In the present writer’s opinion it is, how¬ 
ever, less well suited to ‘Quantum Dynamics’ where predictions concerning measurements 
of events at different space-time points are required. Flavour oscillation experiments are 
of just this latter type. The formulation of QM most naturally adapted to such problems 
is the Feynman Path Integral. 


The story of this approach is really that of the re-introduction of space-time into 
microscopic physics after the quantum revolution of the early 20th Century. Bohr declared 
in his ‘Copenhagen Interpretation’ talk at Lake Como 0 that space and time were 


outmoded classical categories in the quantum world i. The first important step in the 
direction of the rehabilitation of space-time was Dirac’s seminal paper on the Lagrangian 
in QM [|18] , that was the basis Feynman’s formulation of the principles of QM specifically 
in terms of probability amplitudes corresponding to particles moving in space-time 0, |5|. 
The essence of this approach is contained in a single formula, the importance of which, 
for a fundamental understanding of QM, had already been stressed by Heisenberg in 
1929 pfl: 

2 


Pfi = Y,T, 


m l 


H-■J2J2(fm\kn)--{k2\ki)(ki\ii) 


k 2 fci 


(5.1) 


This equation gives the probability of transitions between a group I — Yhi H of initial states 
and a group F = )T m f m of final states. The use of this formula to describe a quantum 
experiment requires knowledge of both the initial (prepared) and final (measured) states. 
Indeed, the corresponding quantum experiment is defined by the experimenter’s choice of 
these states. Only when they are specified is a meaningful comparison of theory with ex¬ 
periment possible. Since typical experimental conditions do not permit to either prepare 
or measure a single quantum state the (incoherent) sums over l and m in (5.1) are neces¬ 
sary. These sums are defined by purely experimental criteria (detector sizes or resolution) 
and are unrelated to microscopic, quantum level, parameters. The states |/c n )...|/c 2 ), \ k\) 
(of which, in the detailed description of any actual space-time quantum experiment, there 
are, actually, a multiple infinity) represent unobserved intermediate states. In the case 
of particle trajectories in space-time, described by space-time coordinates, x tl of interest 
for the present discussion, the phase to be assigned to any given path amplitude, say, 
{fm\xn)■■■{x 2 \xi)(xi\ii) is given by the corresponding classical action divided by h. That 
is, the integral of the classical Lagrangian (divided by h) along the classical particle trajec- 


5 In Reference Jl7]], Bohr wrote: ‘Notwithstanding the difficulties which, hence, are involved in the 
formulation of the quantum theory, it seems, as we shall see, that its essence may be expressed in the 
so-called quantum postulate which attributes to any atomic process an essential discontinuity, or rather 
individuality, completely foreign to the classical theories and symbolised by Planck’s quantum of action. 
This postulate implies a renunciation as regards the causal space-time co-ordination of physical processes.’ 
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tory corresponding to the sequence of space-time points aq, aq ,...x n where t\ < f 2 --- < t n . 
An important feature of Eqn(5.1) is that interference occurs only between amplitudes 
corresponding to different unobserved intermediate states never between different initial 
and final states. 


For the flavour oscillation problem, the general path integral formula (5.1) may be 
simplified considerably. The initial and final states correspond to the space-time posi¬ 
tions i of the source particle and the detection event. As the particles propagate over 
macroscopic distances, in free space, they follow classical (straight line) trajectories. The 
corresponding Feynman path integral then reduces to a Green function, or propagator, 
that is the Fourier transform of the well-known momentum space propagator of a free par¬ 
ticle in quantum field theory. The explicit, general, expression for the invariant space-time 
propagator of a fermion was given in one of Feynman’s early papers on QED |2(| . More 
recently, it has been rederived directly, from the covariant path integral, for an arbitary 
massive particle, by Mohanty [|2T|. Apart from solid angle correction factors, that are not 
relevant to the conventional 1-dimensional discussion of flavour oscillations, the propa¬ 
gator for on-shell particles, or for virtual particles over macroscopic time-like invariant 
intervals such that cAt > 1/mp (where mp is the pole mass of the particle) is given, for 
a stable particl J , up to constant factors, by the simple expression: exp (—imp At). This 
means that it is mathematically identical to the ‘Lorentz invariant plane wave’ of Eqn(2.1). 
However, its physical meaning: a Green function that gives the amplitude for a particle to 
be found at a given space-time point, when, at some prior time, the particle was situated 
at another, well defined, spatial position, is quite different to that of the ‘plane waves’ 
(energy momentum eigenstates) of conventional Schrodinger-Born Wave Mechanics. The 
latter, interpreted as wavefunctions according to the Born prescription, are not square 
integrable, and so are devoid of any information concerning the position of the particle. 
In the case of flavour oscillation experiments, on the other hand, the experimental knowl¬ 
edge that can be obtained concerning the position of the propagating particles, once they 
are created by the decay process, is essentially classical, i.e., unaffected by uncertainties 
of a purely quantum mechanical nature. The meaning of this remark will be clarified by 
a concrete example to be given below. In fact ‘quantum uncertainty’ as manifested in the 
Heisenberg Uncertainty Relations results only in the uncertainty in the time at which the 
propagating particle is produced by the unstable source particle. The associated spread 
in the physical mass of the source particle, and of any unstable recoil particles, then pro¬ 
duces, via energy-momentum conservation in the decay process, momentum or velocity 
smearing of the neutrino mass eigenstates. 


To take a specific example, consider the decay of a pion at rest: tt + —> p + iq. As the 
pion is unstable, it has, in general, a physical mass W n that is different from the most 
likely physical mass, the pole mass, m n . Similarly, the physical mass W fl is different from 
the pole mass m /t . The Breit-Wigner widths, T, of the physical mass distributions are 
related to the particle mean lifetimes, r, by the energy-time Uncertainty Relation Tr = 1. 

6 In the following, quantum states will be specified in ‘space-time’ rather than in ‘configuration space’ as 
is conventionally done. The latter definition is convenient when describing, say, the spatial wavefunction 
of a bound state, such as the hydrogen atom, which requires 6 spatial coordinates for its complete speci¬ 
fication. However, flavour oscillation experiments relate events at different space-time, not configuration 
space, positions. 

7 An unstable particle is described by adding a negative imaginary term to the pole mass: mp —> 
mp — *T/2. 
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Since the initial state of the path amplitudes, that of the source pion, is the same in 
all path amplitudes, the kinematical effects of the non-equality of W n and m n will be 
the same for all path amplitudes, i.e. the resultant velocity smearing of the neutrinos is 
incoherent. On the other hand, since the recoil muon is unobserved, its physical mass may 
be different in the paths corresponding to different mass eigenstates. The corresponding 
Breit-Wigner amplitudes will then create a coherent velocity (or momentum) wave packet 
for each eigenstate. This is a simple application of Eqn(5.1). The spatial trajectory, for 
fixed momentum, p, yielding the path amplitude: 

PA = {f m \x n ,p)...(x 2 ,p\xi,p)(x 1 ,p\ii) (5.2) 

is replaced by the sum of path amplitudes (path integral): 

PI = J2 PA <1 = '52(fm\x n ,p q )...(x 2 ,p q \x 1 ,p q )(x 1 ,p q \i l ) (5.3) 

q q 

As shown in Reference Q, the effect of such wave-packets on the oscillation phenomenon 
is minute, in particular for the case of pion decay, as well as for neutrino oscillation 
experiments in general. As discussed below, it may, however, become more important 
for quark flavour oscillations where the spread in the physical masses of the propagating 
particles is of the same order as their mass difference. 

All the relevant constraints imposed by the Heisenberg Uncertainty Relations have 
been taken into account in the above discussion. In particular, as will be considered in 
more detail in the following Section, there is no ‘spatial wave packet’ associated with the 
neutrinos, and the possibility, (or not) of interference between the path amplitudes does 
not depend upon the existence (or not) of such a hypothetical spatial wave packet. The 
momentum wave packet referred to above is not related, by a Fourier transform, to any 
spatial wave packet with simple physical properties. The muon is an unstable particle with 
an exponential decay law and mean lifetime r M . The Fourier transform of this exponential 
yields, in energy space, a Breit-Wigner amplitude with width parameter T M = l/r M that 
describes the distribution of the physical mass, W^, of the muon. Using energy-momentum 
conservation, the corresponding spread in W M generates a spread in the momentum and 
energy of the neutrino. This spread, weighted by the Breit-Wigner amplitude, generates 
the momentum wave packet referred to previously. Taking the Fourier transform of the 
Breit-Wigner amplitude will give back the original exponential decay law, not a spatial 
wavepacket. 

The irrelevance of the momentum-space Uncertainty Relation to the problem under 
discussion is demonstrated by applying it to the neutrino momentum wave packet in pion 
decay just discussed. With A p u = m^T || the corresponding spatial uncertainty is 
Ax v = 1.27 km. Does this give a limit on the possible knowledge of the spatial position 
of the neutrino? By measuring the time of decay of the pion through detection of the 
decay muon, with the easily obtained precision of 1CU 10 sec, the subsequent position! 
of the neutrino is known with an uncertainty of c x 1CU 10 = 3 cm, since the speed of 
light is precisely known. The known uncertainty in the position of the neutrino is then a 
factor 4 x 10 4 smaller than Ax u as calculated from the Uncertainty Relation. The latter 
evidently has no relevance to the actual or possible knowledge of the spatial position of 

8 It’s direction can be determined by simultaneous measurement of the direction of the decay muon 
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the neutrino. Except for the spread in decay times of the parent pion (described by the 
energy-time Uncertainty Relation) there is then no ‘quantum uncertainty’ on the position 
of the neutrino. Once created, its position can be known with classical precision. 

The final state in Eqn(5.1), in its application to flavour oscillations, is the spatial 
wavefunction of the detection event. At the microscopic level, this can be considered as 
localised to within an uncertainty given by the spatial wavefunction of the target particle, 
which is of atomic dimensions (say 1CT 8 — 1CT' cm). This, however, does not reflect 
the experimental knowledge of the position of the detection event, that is limited, by 
considerations of experimental resolution, to the group of states F in Eqn(5.1). In a 
typical experiment, the experimental resolution might be of the order of a centimetre, 
seven or eight orders of magnitude larger than the microscopic localisation of the target 
particle. Thus, to a very good approximation, the states \f m ) can be assumed to be spatial 
eigenstates. The correction for the actual distribution of detection events is then included 
by performing a suitable average over the oscillation probability, as the contribution of 
different final states is incoherent, requiring classical addition of probabilities. An exactly 
similar argument applies to the initial states \ii) corresponding to the spatial position 
of the source particle at some well defined time. These states also can be assumed, for 
calculational convenience, to be spatial eigenstates. It is important to stress that this does 
not imply that the momentum uncertainty of the source particle is infinite (inconsistent 
with the assumption that it is rest). This is a gross misinterpretation of the meaning of 
the momentum-space Uncertainty Relation. Consider the case that the source particle 
is bound in an atom. The actual quantum positional uncertainty, at the time chosen to 
define the initial state, could be included by writing its spatial wavefunction explicitly in 
the path amplitude. However this distance is so much smaller than the experimenter’s 
actual knowledge of the position of the source, which is all that matters to accurately 
evaluate Pfi, that it is evident that the same result will be obtained if, for simplicity, 
a spatial eigenstate is assumed for the source particle. In an analagous way, energy- 
momentum eigenstates are conventionally assumed for all initial and final state particles 
in the momentum space calculatation of invariant amplitudes in quantum held theory. 
This in no way affects the correctness of the results obtained. A naive application of the 
momentum-space Uncertainty Relation, as done above for the path amplitude calculation, 
to this case would imply an infinite uncertainty on the space-time position of the scattering 
event or decaying particle. Actually the space-time parts of the wavefunctions of the 
incoming and outgoing particles play no role in a momentum space calculation except for 
imposing the constraints of exact energy and momentum conservation. If this is assumed 
ab initio the space-time parts of the wavefunctions can be omitted completely without 
changing the results of the calculations. 

Of course, in any actual experiment, the source particle is not at rest. If it is negatively 
charged and forms an atomic bound state it will have some Fermi motion. If it is positively 
charged, and remains unbound, it will be subject to random thermal motion. Corrections 
for these effects can be included in the calculation of the oscillation frequency by a suitable 
averaging procedure. Since they affect only the definition of the group / of initial states, 
these corrections are incoherent, and, just as in the case of the positions of the source 
particle and detection event, the averaging must be done at the level of probabilities, not 
amplitudes. In fact averaging the source momentum distribution at the amplitude, rather 
than probability, level will result in the generation of a spurious ‘momentum wavepacket’ 
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with parameters characteristic of the initial momentum wavefunction. As discussed in 
the following Section, just such a mistake has been made in many papers on the QM of 
neutrino oscillations in the literature. 

Before defining the intermediate states \kf in Eqn(5.1) for the case of lepton flavour 
oscillations, it is important to consider the nature of the decay process in which the 
neutrino mass eigenstates \vf) and \v 2 ) are produced, as already mentioned in Section 1. In 
the Standard Electroweak Model, these states are created by diagrams containing real or 
virtual W bosons. In the case that neutrinos are massive, and the different eigenstates have 
different masses, there is a non-diagonal, unitary, flavour/mass mixing matrix, the MNS |J 
matrix, that is strictly analgous to the better-known CKM |10| quark flavour/mass mixing 
matrix. The existence of such a non-diagonal matrix implies that lepton (or generation) 
number is not conserved in transitions in the lepton sector just as generation number is 
not conserved in the quark sector because of the non-diagonal nature of the CKM matrix. 
For example, the decays: 

7T+ -> (W + )* -> p + Z/1, VT+ -> (W + )* -> p + U 2 (5.4) 

where (W + )* denotes a virtual W + boson, are strictly analagous to the quark sector 
processes, where an on-shell W + boson decays into quarks: 

W + -»• cs, W + -> cd (5.5) 

in which the first (second) transitions in Eqn(5.5) are Cabibbo allowed (suppressed), 
the neutrino mass/flavour mixing angle 0 in Table 1 plays the same role in (5.4) as 
the Cabibbo angle 6 C in (5.5). Note that the /i + in (5.4) and the c quark in (5.5) are 
both second generation particles, whereas the decay antiquarks are either first or second 
generation. This suggests that the neutrino mass eigenstates should also be associated, 
as are the quark mass eigenstates, to different generations. Just as the W decays in (5.5) 
are independent physical processes, so also are the pion decays in (5.4). This has the 
important consequence, for the following discussion, that |zq) and \u 2 ) may be produced 
at different times in the corresponding path amplitudes. As will be seen, this leads to 
an important new contribution to the oscillation phase. The incoherent nature, as in 
(5.4) above, of decays into different mass eigenstates has been previously pointed out by 
Shrock 0, |. 

What has been assumed until now in all discussions of the QM of neutrino oscillations 
in the literature is that is that lepton number is conserved in pion decay, so that the unique 
physical process is 7r + —> Thus the ‘flavour momentum eigenstate’ [u^) is assumed 

to be produced at a fixed time. It is then rewritten, at this time, in terms of U\) and \u 2 ) 
using the mixing amplitudes given in Table 1. These mass eigenstates subsequently evolve 
in time according to the different propagators discussed in Sections 2, 3 and 4 above. 
The present writer’s opinion is that this procedure is wrong. If neutrinos of different 
generations have different masses, lepton number cannot be conserved at charged current 
weak vertices, just as, in the quark sector, generation number is not conserved at such 
vertices. The non-existence of a consistent theoretical description of ‘flavour momentum 
eigenstates’ has been previously stressed in the literature [^2|. Once the incoherent nature 
of, for example, the two processes in (5.4), is recognised, it will be seen that the conceptual 
problems that have beset the treatment of the QM of neutrino osillations in the past, as 
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evidenced by the many kinematical and geometrical inconsistencies discussed above, all 
disappear. The neutrino oscillation phenomenon occurs even when the constraints of of 
both exact energy-momentum conservation and space-time geometry are imposed. The 
latter is simply the condition that the source particle is defined and the detection event 
observed, each at a unique space time point, and that both the source particle and the 
neutrinos follow classical space-time trajectories that, together, link these two points. 

The detailed application of Eqn(5.1) to the derivation of the oscillation probability for 
a simple two-flavour problem, with the source at rest, will now be discussed. A specific 
example is given by T /t — > z/ e ’i oscillations following the decay, at rest, of a 7r + . The 
initial state corresponds to the undecayed source particle, assumed to be at rest (or to 
have, in any case, a very small and random Fermi or thermal motion), at some arbitary 
time t 0 - At later times t\ or t 2 the source particle may decay into the eigenstates |zq) 
or \v 2 ) respectively. These two possiblities are, of course, in classical physics, mutually 
exclusive. In QM however, as the two different decay modes are indistinguishable, when 
only the detection event is observed, the corresponding amplitudes, not probabilities, must 
be added. The final state corresponds to the detection event, which occurs at time to 
at distance L from the source, and which may be produced by an interaction of either of 
the (indistinguishable) neutrino mass eigenstates with a particle of the detector. Suppose 
now that m\ > m 2 . Since either neutrino mass eigenstate may create the detection event, 
and since v 2 > V\ , it follows that t 2 > t\. 


The different quantum states in Eqn(5.1) may now be identified. The initial state is 
| Si,t 0 ) where S denotes the source particle. The label / corresponds to different source 
positions in the experimental apparatus, as well as, possibly, small random momenta of 
the source particle. The source wavefunction then propagates in time (but not in space) 
over the interval: to < t < U. This evolution is described, in the case of a source particle 
of mass ms and decay width Ts, by the Feynman invariant space-time propagator |20|, pi 


(Si,ti\Si,t 0 ) = exp 


( ims — )(Tj Tp) 


(5.6) 


(source particle of mass ms) 


where r t and r 0 
times ti and t 0 . 
process is more 
propagator is: 


are the proper times in the rest frame of S corresponding to the laboratory 
For the case of a source which is a /3-radioactive nucleus, where the decay 
appropriately described by non-relativistic quantum mechanics, the source 


(Si,ti\Si,t 0 ) = exp 


(~ iE P ~ y)(a " t o) 


(5.7) 


(/? — decay of a nucleus) 


where, neglecting the recoil energy of the daughter nucleus, Ep is the total energy release 
in the /Tdecay process. The formula (5.7) may be derived in the same way as the analagous 
one for a radiative transition in atomic physics, by the use of perturbation theory and the 
time-dependent Schrodinger equation [p3|. The appropriate amplitude to describe the /3- 
decay transition is then the matrix element (Nf\T\N t ) where N t and Nf denote the parent 


9 ' —> ;/ e ’ is only a shorthand notation to indicate that a muon is involved in the production process 

and that a electron is observed in the detection event. This notation is used, without quotes below, as 
is conventional. It should, however, not be forgotten that and v e do not exist as physical particles if 
neutrinos are massive. 
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and daughter nuclei respectively. Since these are in energy eigenstates the Schrodinger 
Equation predicts that they evolve in time as: At,-, t) = At,- , to) exp[—iEj(t — t 0 )], (j = 
i,f). It follows that: 

(N f ,t\T\Ni,t) = (Nf,t 0 \T\Ni,to) exp[-i(Ei - E f )(t - t 0 )] 

~ (N f ,to\T\Ni,t 0 )exp[~iEp(t — to)] (5.8) 

where, in the last member of Eqn(5.8), the recoil energy of the nucleus Nf has been ne¬ 
glected. In the path amplitude formalism, (Nf,t 0 \T\Ni,to) may be considered as a time 
independent amplitude (the same in all interfering path amplitudes) while the exponen¬ 
tial factor in (5.8) represents the propagator of the parent nucleus. It is assumed, for 
simplicity, in (5.8) that t — t 0 -C 1/IV For a source at rest where t = r, Eqn(5.8) leads to 
(5.7). It is interesting to note that, if AT, and Nf in (5.8) are replaced, for pion decay, by 
the quark degrees of freedom, via the correspondences: A r ,) ~ |7r) and \Nf) ~ |0), where 
|0) denotes the vacuum (zero energy) state of the quark sector, then the rest frame pion 
propagator is predicted by Eqn(5.8) to be : exp[— im^Ar], in agreement with the Lorentz 
invariant Feynman space-time propagator of Eqn(5.6) above. 

The source particle then decays into a final state containing zy with the time-independent 
amplitude: 

A(i<-S) = (i\a) (a\T\Si) 

= (i|r|S,> (5.9) 

In the first line of Eqn(5.9), the neutrino mass/flavour mixing amplitude is shown explic¬ 
itly. However, as explained above, like the elements of the CKM matrix, this amplitude 
is an intrinsic part of the charged current coupling, so that the amplitudes A(1 <— S) 
and A( 2 <— S) describe independent physical processes. Thus there is no reason for the 
assumption t\ = f 2 that has been made in essentially all discussions of the QM of neutrino 
oscillations in the literature to date. 


The third element in the path amplitude is the invariant space-time propagator of the 
neutrino ]20], ^1| 

(i, x D , t D \i, Xi, ti) = exp {-imi{r D - r*)] (5.10) 


The fourth and last element in the path amplitude is the time-independent amplitude 
of the detection process: 

A(d^i) = (d m \T\(3)((3\i) 

= ( d m \T\i ) (5.11) 

where the label m specifies the spatial position of the detection event, as well as possibly 
the directions and energies of the particles, in the final state of the event, that are detected. 

The path amplitudes are then given by combining (5.6) or (5.7), (5.9), (5.10) and 
(5.11): 

Ai{(3 <— a) = (dm\T\P){P\i){i,x D ,t D \i,Xi,ti){i\a){a\T\Si){Si,ti\Si,to) (5.12) 

The phase increment A0j can be read off directly from Eqns(5.6) or (5.7) and (5.10): 

A 4>i = 7Tli(T D - T { ) + E S {Ti - T 0 ) (5.13) 
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where Eg = ms or Eg for unstable particle or nuclear 0-decay sources respectively. Now, 


Ti ~ T 0 = 


ti t 0 
IS 


= ti - t 0 + 0(v 2 s ) 


— tjj — to -b 0(vi 


— to — to — L 


1 + — 
. 2 pi 
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(5.14) 


so that, using Eqns(2.5) and (5.14), Eqn(5.13) may be written (neglecting any random 
motion of the target) as: 
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giving for the oscillation phase: 
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(5.16) 


where 

Am 2 = m\ — m\ (5-17) 

It can be seen that, in addition to the contribution due to neutrino propagation, Erst 
given by Gribov and Pontecorvo (the Erst term on the RHS of Eqn(5.16)), there is also 
an important contribution due to the propagation of the source particle during the time 
interval t\ <t < t 2 . For example, in the case of pion decay at rest, Eg/2p 0 = m n /2p 0 = 
2.34, so that the oscillation phase is 34% larger than the Gribov and Pontecorvo result 
and a factor 2.68 times larger than the prediction of the standard formula. In order 
to calculate the oscillation probability, from flavour a to flavour /3, at the fixed source- 
detector separation, L , the probability corresponding to the sum of the amplitudes in 
Eqn(5.12) must be integrated over the detection time to. For this, it is convenient to 
introduce the times-of-flight t { 1 ', tJ/ of the neutrinos: 


t { 1 = t D -U 


1 + — 
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(5.19) 


tmin takes the values t 0 + t{\ to + t { 1 and to + t { 1 for the squared amplitudes for neutrinos of 
masses mi, m 2 and the interference term respectively, where it is assumed that mj > m 2 , 
so that t { 1 > t { 1 The to integral cancels the time-of-flight dependence of the squared 
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amplitude terms and gives a factor exp[—Ts^ —)/2] multiplying the interference term. 
Thus, using Eqn(5.18) the following time-averaged oscillation probability is obtained: 
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(5.20) 


In all conceivable neutrino oscillation experiments the exponential damping of the inter¬ 
ference term is very small. Since observation of the oscillation requires the argument 
of the cosine in Eqn(5.20) to be of order unity, i.e., Am 2 /po ~ 1, the damping term is 
typically ~ exp(— Ys/^Po). In the case pion decay, where Ts = 2.5 x ICE 14 MeV, and 
Po = 29.8 MeV, the damping factor is ~ 1 — 4.1 x ICE 16 . For Am 2 = (leV) 2 , the distance 
L, in pion decay at rest, is 1.9 x 10 16 m or 2.0 light yr for 50 % damping of the interference 
term. 


Two important special cases of Eqn(5.20) are P e/1 , as, for example in —> u e oscilla¬ 
tions following 7r + decay (v e appearence) and P ee as in v e —> v e following /E-decay using 
a nuclear reactor as a source (u e disappearence). Neglecting the exponential damping 
of the interference term, and using the flavour/mass mixing amplitudes in Table 1, the 
following predictions are obtained: 


| (d m I T I e) (fjL I T I Si) 1 2 sin 2 9 cos 2 9 


P ef , = 2 


Am 2 

( m i r ,\ T 

1 — cos- 

-2^-1 \L 

Po 

\ 2 Po ) 


P = 

1 ee 


|<d m |T|e><e|T|S))f 


sin 4 9 + cos 4 9 + 2 sin 2 9 cos 2 9 cos ^ m (—- — 1^ L 

Po V 2 Po J 


(5.21) 

(5.22) 


The derivation of Eqn(5.20) has neglected velocity smearing of the neutrinos due to 
the finite decay widths of the source, and also possibly, the recoil particles. These effects 
have been estimated for the case of pion decay (using a Gaussian approximation for the 
Breit-Wigner amplitudes) in Reference |J. They are found to be more than ten orders of 
magnitude smaller than the already minute correction term for liftime damping discussed 
above. Corrections to Eqn(5.20) due to finite experimental source and detector sizes , 
as well as the effect of thermal motion (for a positive pion source) are also calculated in 
Reference 0. 


An extended discussion of the application of the Feynman path amplitude method to 
flavour oscillations of neutral kaons and b-mesons will be presented elsewhere ]24|]. Here, 
only a few remarks on the salient differences with respect to neutrino oscillations will be 
made. 


(i) Because of the very small fractional mass differences between different neutral kaon 
and b-meson mass eigenstates [p5[| : 
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5.9 x 10~ 14 

(5.24) 
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The Lorentz -7 factors relating the proper time and laboratory times of the mass 
eigenstates are then equal to within the fractional differences quoted in Eqns(5.23) 
and (5.24). Thus, unlike for the case of neutrinos, the ‘equal velocity’ hypothesis 
discussed in Section 2 above, is expected to be a good kinematical approximation. 


(ii) The unobserved propagating particles K s, K l and B^, B h are unstable, with decay 
widths that for K 5 , B^, and B# are of the same order as the mass differences: 

AfflK 

r L 

AfflK 

r B , 

Arris 

(5.28) 


= 2.10 

(5.25) 

= 0.037 

(5.26) 

= f Bh = 1.37 

A777-B 

(5.27) 


Unlike in the case of neutrino oscillations, where the natural widths of the source 
and recoil particles give only very small variations of velocity, the decay widths of 
the K 5 , Bl and B# then produce velocity smearing effects of the same order as the 
velocity differences resulting from the mass differences. Such effects must then be 
taken properly into account in the calculation of the flavour oscillation probability. 
However, because U C T^, the natural width of the may, in first approximation, 
be neglected. Including correctly the kinematical effect of the physical mass, W, 
of the propagating particles modifies the ‘on-shell’ expression (2.5) for the phase 
increment to: 

rriiAti rriiWiAti niiW',1. 

A = rriiATi = -=---=- 

7 i E i Pi 

It is important to note that although, for example, the Kg is significantly off-shell 
on the scale of Atttk, the condition: cAr 1/Ws is still well satisfied in a typical 
flavour oscillation experiment, so the off-shell K 5 can still be considered to propagate 
as a classical particle according to Eqn(5.29). It may also be remarked that the 
pole mass, m*, in Eqn(5.29) results from the Fourier transform of the invariant 
momentum-space propagator, and so is constant for each particle species, unlike 
the variable physical mass W t . 

(iii) Unlike neutrinos which, in terrestial oscillation experiments, are produced inco¬ 
herently as single quantum states from a coherent source (an unstable particle or 
nucleus) neutral kaons may either be produced incoherently, in an incoherent inter¬ 
action, such as: 7 r”p —> A(Kl, Ks) or K~ p —> ii(Kl, Ks), or in correlated ‘entangled’ 
states, from a coherent source, as in: 0 —> K^Ks. To date, B-mesons have usually 
been produced in entangled states from coherent sources, for example Z —» BlBhX 
or T(4S') B l B^. 



Because of the points (ii), (iii) above, the path amplitude analysis of K s — K L and 
E>l — Bh oscillations is more complicated than for the case of neutrino oscillations. The 
standard procedure, to date, in analysing experiments, is to assume equal velocities (that 
seems, in view of Eqns(5.23) and (5.24) to be reasonable) so that, for example, Ary = 
A Ts = At and use the first member of Eqn(5.29), together with the relation: At = 
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rriKoL/p K o to give the formula: 


(5.30) 


A 4>i = 
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Setting tttko = (rriL + ms)/2 in Eqn(5.30) gives, for the oscillation phase: 

Am,K / m K oL (uil — ms)(rriL + rris)L (■ m 2 L — m| )L 
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(5.31) 


This is the same as the standard formula for neutrino oscillations, also derived on the 
basis of the equal velocity hypothesis. At the time of writing, the author is not able to 
comment on the correctness or otherwise, of Eqn(5.31). 


7r 


Taking into account exact energy-momentum conservation in the process: 

"p —> A(K l , K s ), but assuming all particles to be on-shell, the authors of Reference 
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obtained a different formula to Eqn(5.31), as well as predicting correlated spatial oscilla¬ 
tions of in the decays of the neutral kaons and the lambda. 


6 Discussion 


There are two essential differences between the work presented in the present paper 
and References j3|, |BJ, and all previous treatments of the QM of neutrino oscillations in 
the literature: (i) the realisation that the different neutrino mass eigenstates are produced 
in independent physical processes, not as a coherent lepton flavour eigenstate, and (ii), 
the assumption that, in the covariant Feynman path amplitude description, the neutrino 
mass eigenstates follow essentially classical space-time trajectories; there is no spatial 
‘wave packet’ associated with their propagation. Historically, as will be seen, the failure to 
notice (i), leading to the universal assumption that the different neutrino mass eigenstates 
are produced coherently and at the same time, required the introduction of a spatial wave 
packet, since otherwise the neutrino mass eigenstates, produced at the same time and 
propagating with different velocities could not produce the detection event at a well- 
defined space-time point. Only, it was argued, by introducing the spatial ‘fuzziness’ 
associated with the hypothetical spatial wave-packet would interference be possible. In 
the Feynman path description there is no such spatial ‘fuzziness’. The detection event can 
be produced by the interaction of either (indistinguishable) neutrino mass eigenstate with 
the target particle in the detector. Since the space-time trajectories of the neutrinos are 
essentially classical, with different velocities, this implies different production times for the 
neutrinos in the alternative interfering amplitudes. In view of (i) and the long lifetimes of 
typical neutrino sources different decay times are allowed. This is the fundamental reason 
for the possible occurence of the oscillation phenomenon. The propagator of the source 
particle contributes to the path amplitudes associated with each neutrino mass eigenstate, 
and, as a consequence of the different decay times, gives an important contribution to the 
interference phase, as demonstrated by the calculations presented in Section 5 above. 

The important point (i) above was pointed out many years ago by Shrock 0,0 but 
was never, to the author’s best knowledge, applied before the work reported in this paper 
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and References |3|, p] to the neutrino oscillation problem. Indeed, the contrary, incorrect, 
hypothesis of an initial ‘lepton flavour eigenstate’, that is a coherent superposition of 
different mass eigenstates, has, instead, been universally assumed in the treatment of 
neutrino oscillations. 


A brief historical review of the introduction of the ‘wave-packet’ concept into the de¬ 
scription of neutrino oscillations will now be made. The first descriptions of neutrino 
oscillations were based on plane-wave propagators either (tacitly) covariant in the case of 
Gribov and Pontecorvo |]J, or, temporal only, as in the case, for example, of Fritsch and 
Minkowski || who were the first to derive the standard oscillation phase. The wave-packet 
concept in the context of neutrino oscillations was first introduced by Nussinov [^7[], to¬ 
gether with the idea of a ‘coherence length’ for neutrino oscillations. Nussinov discussed 
the effect analagous to collision broadening of atomic spectral lines for neutrino sources 
in the Sun. Both the ‘wave-packet’ and ‘coherence-length’ concepts are related to a clas¬ 
sical wave rather than a quantum mechanical description of the associated phenomena. 
A source of classical waves of finite duration will evidently produce a wave-train of finite 
length. Fourier analysis of this spatial wave-train will result in a spread in the momentum 
of the associated plane-waves described by a coherent momentum wave-packet. This clas¬ 
sical description does not match the sequence of spatio-temporal phenomena underlying 
the line broadening effect as described by quantum mechanics. Actually, the unstable 
source produces particles (photons or neutrinos) in a process that has a characteristic 
timc0 much shorter than the mean decay time of the source. Interference effects will 
occur if some process, intitiated by a decay particle, can correspond to different produc¬ 
tion times of the latter. In the case of, say, atomic line broadening as observed in the 
interference-fringe contrast in a Michelson interferometer, this time difference results from 
the different propagation times of the photon in the different arms of the interferometer. 
In the case of neutrino oscillations, with a source at rest, it results from different times- 
of-flight of the different mass eigenstates over the same spatial distance. In both cases the 
line broadening effect results from perturbation of the coherent source during the time 
interval between the two emission times corresponding to the quantum interference condi¬ 
tion. Thus the fundamental physical parameter governing the interference damping is the 
effective lifetime of the coherent source, due either to its spontaneous decay rate or to col¬ 
lision processes. Nussinov correctly identified this parameter, r e //, as the one controlling 
the damping of the oscillations, but then introduced a hypothetical ‘spatial wave packet’ 
of length CT e ff. This is actually a classical wave analogue of the series of spatio-temporal 
quantum processes just described. As is well-known from atomic radiative transitions 
the classical wave and QM calculation both lead to a Lorentzian line-shape and so are 
equivalent if only a description of the momentum distribution of the produced particles is 
required 0. As shown in the calculations in Section 5, a correct description of interference 
damping requires the spatio-temporal sequence of quantum events to be properly taken 
into account. This is not possible in the classical wave description in which a hypothetical 
spatial wave packet is introduced. It must be stressed that in the quantum mechanical 
description of atomic radiative transitions there is no wave packet associated with pho- 


10 Possible ways to experimentally measure such a ‘p rocess duration’ time are discussed in Reference |28j . 

11 This is only true as regards the mathematical form of the line shape. The width of the distribution 
predicted by the QM calculation depends on the lifetimes of both the initial and final states, so there is, 
in general no direct relation between the line width and the length, cr e //of the analagous classical wave 
packet. 
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ton propagation. The Fourier transform of the exponential decay law of the initial atom 
produces a Lorentzian distribution in the energy of this atom. When it decays by photon 
emission, energy conservation then produces a smearing in the energy (or wavelength) of 
the photon that reflects the energy uncertainties of both the initial and final atoms. But 
the photon propagates in space-time like a particle with velocity c. There is no associated 
spatial, or momentum, wave packet. Indeed, a photon is a particle and not a classical 
wave. 


Krauss and Wilczek [|29| who also estimated the effect of collision broadening on solar 
neutrino sources, although also mentioning the wave packet description, understood that 
neutrino oscillations require different emission times in the interfering amplitudes: 


‘Since these mass eigenstates propagate with different velocities (for fixed energy) the 
desired interference is between neutrinos emitted at different times.’’ 


Actually, the mass eigenstates do not have fixed energy (or momentum) but certainly 
have different velocities, so that the above assertion concerning the interference mechanism 
underlying neutrino oscillations is correct. 


The first extended discussion of the QM of neutrino oscillations in terms of wave pack¬ 
ets was made by Kayser |3(J. In contradiction with the conclusions of Shrock published 


m the previous year [0 the initial state of the neutrinos was assumed to be a superposi¬ 
tion of mass eigenstates with definite lepton flavour, and the propagating neutrinos were 
assumed to have equal momenta and different energies. Thus energy-momentum conser¬ 
vation is violated in the neutrino production process. Kayser claimed, that in order for 
the neutrino oscillation phenomenon to be possible, the uncertainty in the momentum of 
the neutrinos must satisfy the condition: 


A p u Am 2 

A A 9 
Pv Pt 


( 6 . 1 ) 


where A ~ 10 — 100. It was then proposed to realise this condition by introducing a 
hypothetical coherent momentum wave packet with a width consistent with the condition 
(6.1). This condition was derived from the momentum-space Uncertainty Relation as 
follows: It was assumed that the uncertainty, Ax v in the position of the neutrino is 
related to its spread in momentum, Ap U) by the relation: Ax v Ap v > 1. The requirement: 
\Ax v = £ osc , where l osc is the neutrino oscillation length, then leads to Eqn(6.1), it being 
assumed that if Ax u ~ l osc it will be impossible to observe neutrino oscillations. Although 
this latter condition is evidently correct if Ax v refers to the experimental uncertainty in 
the observed position the neutrino detection event, it is easily shown to be completely 
false if, instead, it refers to the theoretical uncertainty calculated using the momentum- 
space Uncertainty Relation. For example, in the case of pion decay at rest the width of 
the coherent momentum wave packet associated with the spread in physical mass of the 
decay muon is A p = = 2.3 x 10” 16 MeV. Taking A = 10, Am 2 = (leV) 2 and 

p v = 29.8 MeV, Eqn(6.1) gives the limit Ap v jp v > 1.1 x 10~ 14 as compared to the width 
of the coherent momentum wave packet in pion decay at rest of Ap u /p u = 7.2 x 10" 18 . 
This is three orders of magnitude lower than Kayser’s lower limit for the possibility of 
neutrino oscillations, (6.1). Yet explicit calculation of the corresponding damping of the 
interference term 0 shows it to be completely negligible. As discussed in Reference 0 the 
momentum spread of the neutrinos due to smearing of the physical mass of the decaying 
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pion, or its thermal motion, are much larger than that associated with the physical mass 
of the decay muon. The two former sources are, however, incoherent, and so have no 
associated wave packets. 

In fact, Kayser assumes that the mass eigenstates are both produced and detected 
at the same times i.e. the equal velocity hypothesis is made, in contradiction with the 
different velocities resulting from the equal momentum and unequal energy hypothesis. 
This procedure is justified by the hope that the effect of this logical inconsistency will be 
annulled by the space-time fuzziness introduced by the hypothetical wave packet: 

‘ The wave packet treatment eliminates the need to make some idealising assumption 
by taking both momentum and energy variations properly into account’ 

The ‘idealising assumption’ referred to is actually the exact conservation of energy and 
momentum that is automatic in all covariant calculations of decay or scattering processes. 
The real purpose of the wave packet is, however, rather to enable the neutrinos, incorrectly 
assumed to be always produced at the same time, but separating spatially due to their 
different velocities, to be detected at the same time. This is clearly impossible when, as 
is in fact the case, they move along essentially classical space-time trajectories. 

Like Nussinov, Kayser introduces a hypothetical spatial wave packet of length cr 
where r is the lifetime of the decaying state. This wave packet, which, as explained 
above, is only a classical wave theory analogue of the finite source lifetime, does not exist 
in the QM calculation. The Fourier transform of the exponential decay amplitude yields 
a Breit-Wigner amplitude describing the distribution of the physical mass of the unstable 
source particle. As discussed in the previous Section, as this mass is a property of the 
source particle, any resulting momentum smearing of the produced neutrinos is then an 
incoherent effect. The only coherent momentum smearing (yielding a momentum, not 
a spatial, wave packet) is that associated with the physical masses of any unobserved, 
unstable, recoil particles. Kayser then argues that the length of the spatial wave packet 
must be much shorter than cr: 

‘If we are interested in a neutrino emitted at time t = 0, but we can learn only that 
the emitter was somewhere in a region of length h , then the amplitude for the emission 
to have ocurred at t = 0 at the various points in this region must be added coherently. 
Thus the neutrino wave packet will have a length d ~ hd 

The false assumption is made that both neutrinos must be created at the same time, 
t = 0. It is then assumed, in contradiction with Eqn(5.1) above, for the case of a neutrino 
source at rest, that the amplitudes of neutrinos created at different spatial positions must 
be added coherently. It corresponds to performing the sums over m and / in Eqn(5.1) in 
a coherent manner, i.e. with J2 m J2i inside the modulus squared. An extended criticism 
of this assumption may be found in Reference 0. In Kayser’s interpretation the spatial 
wave packet is needed to ‘delocalise’ the neutrinos, that, because of the assumption of 
a common production time, become spatially separated due to their different velocities. 
In fact, the neutrinos can arrive at the detector at the same time because they may be 
produced at different times in the alternative histories corresponding to the different path 
amplitudes. There is then no need for the hypothetical spatial wave packet introduced 
by Kayser. There are no physical parameters in the QM calculation governing the size 
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of such a wave packet. It does not exist except in an analagous classical wave theory 
calculation. The space-time structure of particle propagation cannot be described by 
such a calculation. 


Because equal velocities are assumed, i.e. that both neutrinos are produced at one 
space time point and both detected at another, Kayser obtains the standard oscillation 
phase, a result shown to be unmodified |30| by convolution with a narrow momentum 
wave packet. 


Ten years after Kayser’s paper on wave packets in neutrino oscillations the first of 
many papers where spatial wave packets, as suggested by Kayser, were implemented 
using Gaussian functions, was written by Giunti Kim and Lee |12|. 
with the statement: 


This paper starts 


‘If neutrinos are massive particles and mixed, a flavour neutrino is created by a weak- 
interaction process as a coherent superposition of mass eigenstates.’ 

This statement is in contradiction with the findings of Shrock [7j, |8| that, if non¬ 
degenerate massive neutrinos exist, lepton flavour is not conserved by the weak interac¬ 
tion, so that the different mass eigenstates are produced incoherently in different physical 
processes. It implies that, as previously assumed by Kayser (and all subsequent authors 
of papers on the QM of neutrino oscillations) that all neutrinos are produced at the same 
time. Therefore they must satisfy the equal velocity hypothesis of Section 2 above, and 
independently of the introduction, or not, of spatial wave packets, or the assumption, 
or not, of exact energy-momentum conservation in the production process, the standard 
result for the oscillation phase must be obtained. 


The introduction of Reference |T2|] contained a list of issues that it was claimed should 
be addressed in order to provide a ‘complete’ treatment of the QM of neutrino oscillations. 
It is instructive to review this list from the perspective of the work presented in the present 
paper and References [k|, ||. The four issues were: 


‘A necessary condition for neutrino oscillations to occur is that the neutrino source 
and detector are localised within a region much smaller than the oscillation length; 
then the neutrino momentum has at least the corresponding spread given by the 
uncertainty principle 0.’ 


(ii) ‘ The energy-momentum conservation in the process in which the neutrino is created 
implies that different mass-eigenstate components have different momenta as well 


as different energies [14].’ 


(iii) ‘The different mass eigenstates must be produced and detected coherently; this is 
possible only if the other particles associated with the production and decay pro¬ 
cesses have energy momentum spreads larger than the energy momentum differences 
of the mass eigenstates.’ 


(iv) ‘The wave function of the propagating neutrino must be a superposition of the wave- 
functions of the mass eigenstates with proper coefficients given by the amplitudes 
of the processes in which the mass eigenstate neutrinos are produced.’ 
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The following comments (C) may be made made on these points: 


(Ci) The first sentence is trivially correct, but is unrelated to the quantum mechanical 
aspects of the problem. The last phrase, based as it is on the arguments of Kayser 
(Eqn(6.1)) is, as argued above, demonstrably without physical foundation. 

(Cii) This is indeed an essential ingredient of a correct treatment, in QM, of neutrino 
oscillations. A corollary is that the neutrinos, propagating over macroscopic dis¬ 
tances, do so as classical particles [^1|]. This condition is violated by the universal 
equal velocity assumption of Eqn(2.12). 

(Ciii) Properly interpreted, the first sentence is correct. It means that the amplitudes 
describing the temporal sequence of events: a) propagation of the source, b) decay 
of the source, c) propagation of a neutrino mass eigenstate and d) production of 
the detection event must be added coherently. It does not mean that the mass 
eigenstates are part of a coherent ‘flavour eigenstate’ in either the detection or 
production processes. The last phrase, as shown explictly for the case of pion decay 
at rest in Section 5 above, is false. 


(Civ) The ‘wavefunction’ referred to here is presumably a lepton flavour eigenstate prop¬ 
agating in space-time. As previously pointed out by Shrock il and shown later 
by Giunti, Kim and Lee themselves ^2|, such a state does not exist. 


The effect of the Gaussian wave packet introduced in Reference ]TIJ on the oscillation 
phase is described in Section 3 above. For a critical discussion of the ‘spatial coherence 
length’ and ‘momentum damping factor’ generated by such wave packets see Reference 
The point of view, explained there, of the present author is that, as the coherent spatial 
wave packet from which they are derived does not exist, they also are without any physical 
foundation. The oscillations are damped by the source lifetime, as shown in Eqn(5.20), 
and also by the width of the coherent momentum wave packet associated with the dif¬ 
ferent possible physical masses of unobserved recoil particles produced in the production 
process. As shown in Reference [[|, both effects are expected to be vanishingly small in 
any forseeable neutrino oscillation experiment. 

An influential paper on the QM of heavy quark flavour oscillations was written by 
Similar ideas were applied to neutrino oscillations in Reference P3 and in an 

0 


Lipkin 
unpublished pre-print 


16]. The starting point of Reference 


was the correct observa¬ 


tion that all experiments measuring flavour oscillations actually observe only a spatially 
varying interference effect. This implies that all decay and detection times should be 
integrated over in order to derive the quantum mechanical probability to be compared 
with experiment. This is done, for example, in the derivation of Eqn(5.20) above. Lipkin 
interpreted this correct statement about the nature of the experiments as implying that 
time should not appear at all in the theoretical description of the experiments. Formulae 
containing the time were referred to as describing ‘non experiments’. In all three papers 
cited above the initial state is incorrectly assumed to be in a pure flavour eigenstate that 
is a superposition of mass eigenstates. Thus, for example the Ks and the Kl, or the 
different neutrino mass eigenstates are assumed to be produced at the same time. Since 
they must be detected at the same time, as only one detection event is observed, the equal 


















velocity hypothesis (2.12) is thus assumed. Since the additional assumption of equal en¬ 
ergies and different momenta of the particles is made (of course in logical contradiction 
with the equal velocity hypothesis) the temporal part of the Lorentz invariant plane wave 
does not contribute to the oscillation phase. The latter is then entirely determined by 
the spatial part, as described in Section 4 above. The standard result for the oscillation 
phase is then obtained. The arguments given by Lipkin to justify the choice of equal 
energies and different momenta of the propagating mass eigenstates are unconvincing. In 
the process K~p —> K n Lipkin states that ‘Energy conservation requires the K to have 
a definite energy. When it is split into Kl and Kg components with different masses the 
two states have the same energy but different momenta’. There would seem to be no 
physical justification for this apodictic statement. Momentum conservation requires the 
K to be produced with a definite momentum. Why should it not then ‘split.’ into two 
states with the same momentum, and different energies? Later Lipkin gives an argument 
based on non-relativistic kinematics (Eqn(2b) of Reference [[OJ]) as applied to the pro¬ 
cess K~p —> K°n to justify the neglect of the energy difference between the Kl and Kg. 
Repeating the calculation using relativistic kinematics, as more appropriate to typical 
experimental conditions, shows instead that A p/p = AE/E, so that there is, in this case, 
no kinematical justification for the equal energy hypothesis. Actually the Kl and Kg, like 
the different neutrino mass eigenstates, are produced incoherently, in different physical 
processes, so that there is no production of the state ‘K ’. Conservation of energy and 
momentum then shows that they must have different energies, different momenta and 
different velocities. 


In Reference JTl| a calculation of the oscillation phase for the B°B° system is performed 
in the laboratory system using temporal evolution, and allowing different propagation 
times for the different mass eigenstates (Eqn 9b of Reference [[□]]). As discussed in Section 


4 above, this corresponds to the full Lorentz invariant phase in the non-relativistic limit 
where pi <C m*. In this limit the complete 0(m 2 ) Gribov-Pontecorvo result, with the 
oscillation phase a factor two larger than the standard result, is obtained. Lipkin noticed 
this difference but rejected the correct result given, in the appropriate kinematical limit, 
by his Eqn 9b on the grounds that, as the time appeared explicitly in its derivation, it 
corresponded to a ‘non-experiment’. 


In Reference [15], similar arguments are applied to the neutrino oscillation case. Again, 


equal energies and different momenta and the (contradictory) equal velocity hypotheses 
are assumed, leading to the standard oscillation phase. The initial state is required to be 
a superposition of different mass eigenstates with pure flavour. For example, in the case 
of say 7r + —> p + v, the probability of detecting a v e is zero, and of detecting a is unity, 
at L — 0. These conditions are used to fix the coefficients of the mass eigenstate super¬ 
position at L = 0. However, as is shown by inspection of Eqn(5.21) above, exactly the 
same boundary condition is respected by the result of the Feynman path amplitude cal¬ 
culation where the n + decays incoherently into the different mass eigenstates at different 
times and no unphysical ‘flavour eigenstate wavefunction’ is introduced. It was claimed 
in Reference [[L5[ that the energy-momentum and space-time descriptions of flavour oscil¬ 
lations are ‘complementary’ and that including them both leads to ‘double counting’ of 
the oscillation phase by a factor of two. In fact, the Lorentz invariant oscillation phase 
contains (except in the rest frame of the propagating particle) both energy-momentum 
and space-time contributions that must both be included to obtain the correct result. 
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The detailed considerations of Sections 2-4 above show that attempts to use only energy- 
momentum or space-time descriptions in inappropriate kinematical regions leads instead 
to ‘half counting’ of the correct Lorentz invariant phase. 


In Reference [16], Lipkin justified the ‘equal energy’ hypothesis by reference to a paper 


by Stodolsky |3l| which attempted a non relativistic density matrix description of flavour 
oscillations. It was proposed to use ‘stationary’ beams of fixed energy to describe the sys¬ 
tem of propagating mass eigenstates. In this way the introduction of spatial wave packets 
was avoided. The present writer’s opinion is that such an approach is ill-founded. In 
fact, the propagating particles describe classical trajectories, the detailed spatio-temporal 
structure of which is essential for the correct QM description of the phenomenon. This 
information is not available in the non-relativistic density matrix approach, which in any 
case, is not appropriate to describe ultra-relativistic neutrinos. 

The ‘stationary’ source and target description with equal energies and different mo¬ 
menta for the neutrinos as well as Gaussian spatial wave packets for both the source and 
detector was used more recently by Ioannisian and Pilaftsis |32|. As the equal velocity 


assumption (2.12) was also made the standard oscillation phase was obtained. 


It is interesting to note that, in an earlier paper, written together with Kayser 


Stodolsky proposed a covariant Feynman path amplitude approach, akin to that of the 
present paper and References |3], [|, to the description of ‘entangled’ systems such as 
0 —> KsKl. In this case, in contradiction to Reference [[H|], and as previously assumed 
by Kayser H. equal momenta and different energies were proposed. The treatment 
of Reference |33|] differs from that of the present paper and References ||, |6|] in that the 


equal velocity hypothesis (2.12) was assumed, so that the standard oscillation formula was 
obtained and the contribution to the oscillation phase of the coherent source was neglected. 
Also, unphysical spatial wave packets were introduced and the potentially important 
velocity smearing effects due to variations of order A ijik in the physical mass of the Ks, 
as discussed in Section 5 above, were neglected. The present writer is in agreement with 
the main conclusion of this paper, that the ‘collapse of the wavefunction’ often discussed 
in connection with entangled states, is only the collapse of a mathematical abstraction, 
that is actually irrelevant to the QM description of the experiment. It is perhaps, however, 
going a little too far to state, as in the last sentence of Reference [^3| that: ‘The best 
answer, finally to the ‘question of the collapse of the wavefunction’ is that there is no 
wavefunction.’ Indeed, the QM of flavour oscillations is better described in terms of 
Feynman path amplitudes, as proposed in Reference |33|, However, the wavefunction 
does remain an important and useful concept in the description of ‘static’ bound systems 
such as the hydrogen atom. 


Following Rich [ 341 several authors 


have used a non-relativistic Wigner- 
Weisskopf type formalism to describe the complete production-propagation-detection pro¬ 
cess. The propagating virtual neutrinos are assigned the same energy and different mo¬ 
menta. In all cases the equal velocity assumption is made leading to the standard os¬ 
cillation phase. The present writer’s opinion is that such treatments take properly into 
account neither the ultra-relativistic nature of the propagating neutrinos nor the sequence 
of spatio-temporal production and detection events necessary for a correct calculation of 
the oscillation phase. In one paper using this non-relativistic approach [[H]], the effect of 
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the source lifetime was discussed. Although it was correctly concluded that the damping 
effect due to the finite source lifetime is negligible, it was proposed that the characterstic 
momentum spread in the Gaussian wave packet, related to the source, for muon decay 
at rest (actually the Fourier transform of a Gaussian spatial wave packet as proposed by 
Giunti Kim and Lee ]12| following the suggestion of Kayser |I0|) should be ~ 10~ 3 MeV. 
This is the typical momentum due to thermal motion at room temperature. As discussed 
above, any momentum smearing due to this source or to the physical mass of the decaying 
muon is incoherent and not associated with any wave packet. Indeed, in the case of muon 
decay, since all recoil particles are stable, then, unlike in the case of pion decay, there is 
is no momentum wavepacket associated with neutrino propagation. 

The equal energy, different momentum, hypothesis was also made by the authors of 
Reference This was justified by assuming that both production and detection of 

neutrinos resulted from inelastic scattering on an infinitely heavy target. Although such 
an assumption guarantees the kinematical correctness of the equal energy hypothesis, it 
evidently does not correspond to actual neutrino experiments where the neutrinos are 
produced by the decay of an unstable source, and where energy-momentum conservation 
always requires (see Eqns(2.6),(2.7) above) that both momenta and energies and, hence 
the velocities, are different. The equal velocity hypothesis (in contradiction with the 
different momenta of the neutrinos) was also made so that the neutrino propagator is 
purely spatial, giving, as shown in Section 4 above, the standard oscillation phase. 

Many of the features of the covariant path amplitude calculation of Section 5 above 
and References I 1 have been previously introduced into the discussion of neutrino 
oscillations. For example, the Lorentz invariant Feynman propagator for the neutrinos 

imrnM- The only author to introduce explicitly 


has been used in References 21, 32 


the invariant propagator of the source particle was Campagne |I9| . As the equal velocity 
hypothesis was also made there is no contribution to the oscillation phase from this 
propagator and the standard result was obtained for the oscillation phase. 


The authors of Reference [|43|] recognised that the derivation of the standard formula 
required the equal velocity hypothesis, and that if exact energy-momentum conservation 
is imposed, so that the neutrinos have different times-of-flight, an oscillation phase a 
factor of two larger is obtained. In spite of noting that neutrinos of widely differing 
masses, as expected theoretically, cannot have equal velocities, the use of the equal velocity 
hypothesis was, nevertheless, recommended. A short note of Okun and Tsukerman [fh| 


pointed out the kinematical impossibility of the equal velocity assumption for neutrinos 
with different masses. 


A paper by Giunti |45| considering the analogy between the interference effects in the 
Young double slit experiment and in neutrino oscillations, claimed to demonstrate that 
the extra factor of two in the oscillation phase, due to neutrino propagation, obtained 
when the different neutrino velocities are correctly taken into account & [43|, is incor¬ 


rect. This argument was based on the correct observation that in the Young double slit 
experiment with photons, photon propagation gives no contribution to the interference 
phase. That this must be so is evident by setting m, = 0 in Eqn(5.9) above. Giunti 
then claimed to have rejected the different velocity hypothesis [£I3| and the path ampli¬ 
tude calculations of Reference [H by reductio ad absurdum since, in the analagous Young 
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double slit experiment, a vanishing interference phase, clearly excluded by experiment, is 
apparently predicted. However, Giunti neglects the contribution to the interference phase 
of the coherent source (the excited atom that produces the photon). The contribution of 
the source to the interference phase is, using the formula analagous to Eqn(5.7) above, 
for an atomic radiative transition: 




= E’O, - f 2 ) ~ - i 2 ) = Ey(r 2 - r t ) 


( 6 , 2 ) 


where E*, and E 1 are respectively the atomic excitation energy and the photon energy, 

This formula is identical to Eqn(8) of 


while fi, f 2 and ri,r 2 are defined in Reference [[45 
Reference 


, and since the contributions from the photon propagators vanish, the path 
amplitude calculation gives the usual result obtained in the classical wave theory of light. 
Having wrongly concluded that the photon is produced at a unique time, Giunti introduces 
a hypothetical Gaussian spatial wave packet (unphysical, because it has no relation to the 
physics of the photon production process) and demonstrates that the space-time smearing 
provided by the wave packet allows to recover the same interference phase as in Eqn(6.2) 
above. Also the interference term is found to be damped by a factor dependent on the 
length of the wave packet. As discussed in Reference [[J, as well as above in the present 
paper, this damping factor, derived from the spurious spatial wave packet is also without 
any physical foundation. The damping of the interference term is actually produced by 
the hnite lifetime of the coherent source that limits the value of t\ — t -2 in Eqn(6.2). 

A recent paper by De Leo, Nishi and Rotelli [£I(| has also considered the effect on 
the oscillation phase of different kinematical assumptions. As in Reference El it was 


realised that only in the case of the equal velocity hypothesis is the standard oscillation 
phase obtained; in any other case the phase is a factor of two larger. This agrees with 
the conclusions of Sections 2, 3 and 4 of the present paper. The difference between the 
kinematical discussions of the present paper and those of Reference |I(| is that in the 
former case, in accordance with the experimental conditions, and as previously pointed 

m 


out in Reference 
detection event. 


a constant distance, L, is assumed between the source and the 


In the latter case this distance is allowed to vary, and also different 
creation times are allowed for the mass eigenstates. In the general kinematical discussion 

H 


of Section II of Reference 


which may be compared to Section 2 above (since both 
discuss the Lorentz invariant oscillation phase) it is assumed that At is the same, but 
that Ax is different for the different neutrinos. This is in complete disagreement with the 
experimental conditions of typical neutrino oscillation experiments. In order to permit 
different source-detector separations, a Gaussian spatial wave packet was introduced. It 
was concluded that by a suitable choice of creation times a pure flavour eigenstate can 
be obtained at creation. Indeed the necessary existence of such a state is the initial 
hypothesis on which the arguments given in the paper are based. It may be commented 
that, firstly, as pointed out long ago by Shrock i! no such flavour eigenstate exists since 
the neutrinos are created in separate, incoherent, processes, and secondly, the distance 
L must be constant, since the source is assumed to be at rest. Finally, although it 
is true that the different neutrinos may be created at different times, the spatial wave 
packet introduced to allow the possibility of different source-detector distances does not, 
as argued above, have any physical basis. 


A very recent paper by Beuthe |42| makes the same basic assumptions as an earlier 
paper by Giunti et al j|7j . Both equal energy and equal energy hypotheses are considered, 
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but the equal velocity condition (2.12) is always assumed so that the standard oscillation 
phase is always recovered, and the important contribution to the oscillation phase of 
the coherent source particle is neglected. There is a lengthy discussion of the effects 
of hypothetical Gaussian wave packets used to describe both the source and detector 
particles. 


Further critical discussion of different treatments in the literature of the QM of neu¬ 
trino oscillations may be found in Reference 0. The same paper also describes briefly 
two atomic physics experiments, the ‘quantum beat’ experiment |48|] and the ‘photode¬ 
tachment microscope’ [49, 50, 51] where the Feynman path amplitude description has 


been successfully tested in experiments where spatially varying interference effects, very 
similar to particle flavour oscillations, have been observed. 


7 Summary and Outlook 


The kinematical and geometrical discussion of Sections 2, 3 and 4 above shows that 
the standard formula (1.5), for the contribution of neutrino propagation to the oscilla¬ 
tion phase, is a consequence of the equal velocity hypothesis where it is assumed that 
both interfering neutrinos are produced at the same space-time point. This hypothesis 
is incompatible with the propagation of the neutrinos along classical space-time trajec¬ 
tories if they have different masses. The oscillation phase, calculated at 0(m 2 ), using 
the equal velocity hypothesis is found to be a factor of two smaller than the result first 
obtained by Gribov and Pontecorvo [jl[]. The latter is obtained by imposing both exact 
energy-momentum conservation and a consistent geometrical propagation in space-time, 
taking properly into account the different neutrino velocities. It is referred to above as 
the ‘exact’ 0(m 2 ) formula. Thus the standard formula neglects numerically important 
0(m 2 ) terms as compared to the exact one. 

This conclusion remains the same whether the neutrino propagator is described in 
a Lorentz invariant manner using plane waves (Section 2), whether convolution with a 
Gaussian wave packet is performed (Section 3) or whether only temporal or spatial evolu¬ 
tion of the neutrino wave function is considered (Section 4). In contrast, the assumptions 
of equal momenta and different energies or of equal energies but different momenta give 
only negligible 0(m 4 ) corrections to the exact formula, or to the standard formula, as 
compared to the phase calculated assuming exact energy-momentum conservation. 

The Feynman path amplitude calculation of Section 5 is based on the use of the exact 
Lorentz invariant plane wave formula (or, equivalently, the invariant Feynman space-time 
propagator) for the contribution of neutrino propagation to the oscillation phase. Since the 
different neutrino mass eigenstates are produced in independent physical processes SI-1 
the decay of the source can occur at different times in the amplitudes describing the 
propagation of different neutrinos. The oscillation phenomenon then occurs respecting the 
constraints of both exact energy-momentum conservation and exact space-time geometry. 
The different decay times of the source in the interfering amplitudes then lead to an 
important contribution to the oscillation phase from the space-time propagator of the 
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source. 


Damping of the interference term due to the finite source lifetime or momentum smear¬ 
ing related to the off-shell nature of the source or recoil particles has been previously found 
to be completely negligible in all foreseeable neutrino oscillation experiments ||. 

All previous calculations in the literature in which the standard oscillation phase is 
obtained assume the production of a coherent ‘neutrino flavour eigenstate’ in the decay 
process, that is a quantum superposition of mass eigenstates. This enforces equal pro¬ 
duction times for the all mass eigenstates and hence the equal velocity condition. As 
pointed out in References Bum! such a ‘neutrino flavour eigenstate’ does not exist. 
The different neutrino mass eigenstates are produced incoherently in different physical 
processes. This universal, incorrect, assumption concerning the nature of the initial state 
of the neutrinos thus explains why the standard oscillation phase has been, hitherto, uni¬ 
versally obtained. Thus the contribution to the oscillation phase of neutrino propagation 
has been generally underestimated by a factor of two, and the important contribution 
to the phase of the coherent source (resulting from different source decay times in the 
interfering amplitudes) has been universally neglected. 


As discussed in some detail in Section 6, in order to enable the interference phe¬ 
nomenon, leading to ‘neutrino oscillations’ to occur when the different mass eigenstates 
are produced at the same time inside the ‘neutrino flavour eigenstate’ it was proposed [3Q 
to introduce a spatial wave packet to delocalise the neutrinos. Moving with different ve¬ 
locities along classical trajectories, and produced at the same space time-point, it is clear 
the neutrinos can never arrive together at the unique space-time point of the detection 
event. The spatial ‘fuzziness’ introduced by the hypothetical wave packets, associated 
with each neutrino, was conjectured to enable them both to have non-vanishing ampli¬ 
tudes at the position of the detection event, so that neutrino oscillations can occur. In 
fact, because the interfering neutrinos can be produced at different times there is no need, 
to produce neutrino oscillations, for the spatial fuzziness introduced by the hypothetical 
wave packets. 


The hypothetical wave packet suggested in Reference |10| does not exist in the QM 
calculation of the neutrino production process. The neutrinos (particles) are produced in 
space-time according to an exponential decay law. A ‘wave packet’ is only a (very loose) 
analogy in classical wave theory to the effect of the lifetime of an unstable particle on the 
energy (or mass) distribution of its decay products. In the analogy this distribution is 
given by the Fourier transform of the spatial wave packet, but no there is no such wave 
packet in the QM calculation itself. Thus the (mathematically convenient) Gaussian wave 
packets that abound in the literature on the QM of neutrino oscillations have no physical 
foundation within QM. The relevant (related) parameter, in the quantum mechanical cal¬ 
culation, is the mean lifetime of the source. There is no physical connection between this 
parameter and the length of a loosely analagous wave packet with an arbitary Gaussian 
form. The introduction of such wave packets in the QM calculation mixes up in a confused 
way concepts from QM and a classical wave theory that contains no information on the 
space-time evolution of particle positions. 


The most important additional features of the present paper as compared to Ref¬ 
erences [[|. 6] are, first, the realisation that the standard oscillation phase follows only 
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from the manifestly unphysical equal velocity hypothsis and is quite unrelated to the 
use of Gaussian wave packets or any of the other kinematical assumptions made in the 
derivations. The misleading impression may have been given in References [[3[ || that the 
standard oscillation phase, when obtained in covariant calculations, was a consequence of 
the use of wave packets. The second is the realisation that the incoherent nature of the 
neutrino production process, which is the physical basis of the calculations presented in 
References SB and Section 5 above, was aleady pointed out more than twenty years 
ago by Shrock 0,1 in the published literature. I was not aware of this work at the time 
of writing References [[| SJ. 


In closing a few remarks are made on phenomenology and experimental tests. The 
mass difference, Am^ tand , derived from experimental results using the standard formula, is 
related to that, A m 2 FPl given by the Feynman path amplitude calculation, by the formula 
(valid for any source at rest): 

AmJ. p = CAmL, = (7.1) 


while for neutrino oscillations following two body decays in flight of nltra-relativistic 
charged pions or kaons [[|: 


a _2 _ ( 1: «' R m) Am2 stand 

laUI fp — 

Z 1 l m 


(7.2) 


where R m is defined after Eqn(2.8). In the case of experiments involving neutrino pro¬ 
duction in pion, kaon or muon decay the conversion is straightforward. However, for 
neutrinos produced in nuclear reactors the effective oscillation phase will require a suit¬ 
able average, with appropriate weighting factors, over the decays of all /3-unstable nuclei 
contributing to the neutrino flux. For a given nuclear species the correction factor, C, 
to the standard oscillation phase is given, by Eqn(7.1), as \/{Egjp v — 2), where Eg is 
the total energy release in the decay. This is evidently an immense undertaking for any 
actual reactor-based experiment. Any phenomeological conclusions hitherto drawn from 
the results of such experiments, using the standard formula, must therefore be discarded 
if the oscillation phase is correctly given by the path amplitude calculation. 


Of particular interest, in view of the recent results of the Ka.miokande 
SNO S collaborations, are the /3-decay processes: 8 B 


*Be + e H 


^ and 
+ (zq, zq, zq) con¬ 


tributing to the flux of high energy solar neutrinos. The correction factor, C , in Eqn(7.1), 
is plotted in the range: 0 < p v /E s < 0.45 in Fig.l. For 8 B /3-decay, E s — 28 MeV. 
The correction factor is unity when p v /Eg = 1/3, corresponding to p u ~ 9.3 MeV. 
Near the kinematical end-point at p v ~ 14 MeV, C rises steeply, reaching a maximum 
value of about 1500. Thus neutrino oscillations, observable when C ~ 1, are strongly 
suppressed@ in the near end-point region. For the case of the electron capture reactions: 


7 


Be 


'Li + (zq, zq, zq) contributing line spectra to the solar neutrino flux, p u ~ E* 


where E* is the excitation energy of the unstable 'Be atom and C = — 1, so that the oscil¬ 
lation phase of the path amplitude calculation is the same as that given by the standard 
formula. 


12 i.e. they correspond to a vanishing mass difference or an infinite wavelength for the oscillation in the 
standard formula. 
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As already mentioned in Reference @ , evidence for neutrino oscillations in short base¬ 
line experiments such as LNSD |54|] and KARMEN |[5) can be confirmed or invalidated 
by a search for muon oscillations following pion decay at rest, since, the muon oscillation 
phase is found || to be identical to that of the associated neutrinos given in Eqn(5.21) 
above. Evidently the event rate in such muon oscillation experiments can exceed that 
possible in the search for the associated neutrino oscillations by many orders of magni¬ 
tude. 
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